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Abstract 



The structures of free-standing zirconium nanowires with 0.6—2.8 nm in 
diameter are systematicahy studied by using genetic algorithm simulations 
with a tight-binding many body potential. Several multi-shell growth se- 
quences with cylindrical structures are obtained. These multi-shell structures 
are composed of coaxial atomic shells with the three- and four-strands heli- 
cal, centered pentagonal and hexagonal, and parallel double-chain-core curved 
surface epitaxy. Under the same growth sequence, the numbers of atomic 
strands in inner- and outer-shell show even-odd coupling and usually differ 
by five. The size and structure dependence of angular correlation functions 

and vibrational properties of zirconium nanowire are also discussed. 
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In the past decade, there has been tremendous interests on ultrathin metal nanowires 
from both fundamental low- dimensional physics and technological applications such as nano- 
electronic devices0^i^. Most of the previous studies are based on the tip-surface contact0~i or 
mechanically controllable break junctionsiS, which can be seen as very short metal nanowire. 
In recent experiments, Takayanagi's group has successfully fabricated stable gold wires with 
various diameter of sufficient length suspended between two standsi^i2l, and the helical multi- 
shell structures are observed in those ultrathin wire On the theoretical side, the noncrys- 
talline structures, melting behavior, and electronic properties of ultrathin free standing Pb, 
Al, Ag nanowires have been investigated by Tosatti and co-workers0Hll. By using molecular 
dynamics-based (MD) genetic algorithm (GA) simulations, our group has studied structural, 
vibrational, electronic, magnetic properties of gold, titanium, and rhodium nanowiresEHll. 
However, the current knowledge on the detailed structural characters and growth sequences 
of the metal nanowires is still quite limited. In this paper, we report systematical multi-shell 
structural growth sequence of zirconium nanowires and their vibrational properties. 

In our simulations, the zirconium nanowires with diameters from 0.6 to 2.8 nm are 
modeled by a supercell with one-dimensional (ID) periodical boundary condition along the 
wire axis direction. For the most zirconium nanowires studied, the length of 1-D supercell is 
chosen to be 1.413 nm, which is a reasonable compromise between discovering the helicity in 
ID direction and avoiding the nanowires breaking into clusters upon relaxation. To match 
the periodicity of the multi-shell structures with pentagonal symmetry, several different 
supercell lengths, i.e., 1.177 nm for S5-1, S5-2 wires (see Fig.l), 1.648 nm for S5-3, S5-4 wires, 
have been used. The "effective diameter" , which can be controlled by adjusting the number 
of atoms in the supercell, is used to denote the cross section area of the nanowiresEl'llHl^. 
To search the most stable structure of nanowire, we adopt the genetic algorithm (GA) 
based on MD relaxation^Hll. The vibrational densities of states can be then calculated by 
diagonalizing the dynamical matrix for the optimized configurations. 

The interaction between zirconium atoms is described by a tight binding (TB) many body 
potentiaiii, which has successfully reproduced the high-temperature hcp-bcc transition i 
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bulk zirconiumll^ll^. To further check the vahdity of TB potential in low-dimensional systems 
with reduced coordination number (CN), we have performed density functional calculations 
on small zirconium clusters (Zry, Zria, Zris, Zrig) and bulk zirconium solid by using a DMol 
package^ within local density approximation (LDA). In both TB and LDA calculations, 
the cluster structures are optimized and the similar equilibrium configurations are obtained. 
Typically, the discrepancy for interatomic distance obtained from tight-binding potential and 
LDA calculation are within 0.1 A. The comparison of LDA and TB results on the binding 
energies of small clusters with equilibrium structures are given in Table I. The satisfactory 
agreement between LDA and TB results on Zr clusters in Table I demonstrates the validity 
of the TB potential in the other low- dimensional systems like nanowire. 
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TABLES 



Table I. Comparison of binding energies from tight-binding (TB) potential and density 
functional calculation with local density approximation (LDA) for the small Zr„ clusters 
(n = 7, 13, 15, 19). The binding energies of clusters are given by the percentage of bulk 
cohesive energy and the average coordination number (CN) are presented. The TB results 
on larger Zrag and Zrss clusters are also included in the Table. The equilibrium configurations 
for the Zry, Zria, Zri5, Zrig, Zrag, Zr55 are pentagonal bipyramid, icosahedron, hexagonal 
anti-prism, double icosahedron, truncated octahedron, icosahedron respectively. 
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FIGURES 




FIG. 1. Morphologies of structural growth sequences in zirconium nanowires 
right part for S4-x,S6-x and left one for S3-x, S5-x, and S2-x (x=l-4). 



From the GA simulation, we have obtained numerous isomers for the zirconium 
nanowires. Among these structural candidates, we found that the ones with closed-shell 
structures are usually the lowest-energy configurations at a given wire diameter. The energy 
difference between the lowest-energy closed-shell structures and other open-shell structures 
is about 0.005eV/atom. For example, the binding energy of helical closed-shell S6-4 wire 
(see Fig.l) is -5.829 eV/atom, while a low-energy non-helical open-shell configuration at the 
same size has a binding energy -5.824 eV/atom. In our simulations, we have tested differ- 
ent supercell length and found that the multi-shell arrangements and the helicity of magic 
closed-shell nanowires do not sensitively depend on the chosen supercell length. On the other 
hand, those non-magic open-shell wires are either crystalline- like or irregular non-crystalline 
and have higher binding energies. These results indicate the relative stability of magic 
closed-shell structure and the possibility of observing these metal nanowires experimentally. 

Shown in Fig.l are five kinds of growth sequences of zirconium nanowires obtained from 
GA optimization. These multi-shell nanowires are composed of coaxial cylindrical atomic 
shells. Each shell is usually formed by a number of atomic strands winding up helically 
side- by-side. The pitches of the helix for the outer and inner shells are different. The lateral 
surface of the cylindrical shell forms a triangular atomic network. Such helical multi-shell 
structures have been theoretically predicted for Al, PlU, Aulli, TiS, and Rh nanowiresS, 
and experimentally observed in Au nanowire In Fig.l, we introduce the index S3-x, 
S4-X, S5-X, S6-X, and S2-x (x=l-4) to characterize the multi-shell (x: number of shells) 
growth patterns based on trigonal (S3-x), tetragonal (S4-x), centered pentagonal (S5-x), 
centered hexagonal (S6-x), and double-chain-core (S2-x) packing, respectively. These twenty 
structures shown in Fig.l represent five kind of continuous structural growth sequences 
of zirconium nanowires. Such systematic multi-shell structural growth patterns in metal 
nanowires have never been reported before. 

With the notation of n-nl-n2-n3-n4, we can describe the closed-shell magic nanowires 
consisting of coaxial shells, each of which is composed of n, nl, n2, n3, and n4 helical atomic 
strands from outer to inner (n > nl > n2 > n3 > 'ri4)0S0. As shown in Fig.l, 



wires 



S3-l(3), S3-2(8-3), 83-3(13-8-3), and 83-4(18-13-8-3) constitute the sequence of trigonal 
multi-shell packing. The structure of the thinnest wire 83-1 (diameter D=0.574 nm) is 
formed by three-strands packing helically (see Fig.l). 8imilarly, wires 84-1(4), 84-2(9-4), 
84-3(14-9-4), and 84-4(19-14-9-4) constitute growth patterns with one-, two-, three-, and 
four-shell tetragonal packing. This sequence starts from the single-shell wire 84-1 with four- 
strands helical packing. In both the trigonal and tetragonal multi-shell growth sequence, 
the outer and inner shells differ by five helical atomic strands. 

In contrast to the helical structures found in 83-x, 84-x wires, the wires 85-1(5-1), 85- 
2(10-5-1), and 85-3(15-10-5-1) under the centered pentagonal multi-shell sequence are not 
helical, while 85-4(20-15-10-5-1) come back to helical. These structures can be viewed as 
deformed icosahedral packing. Wires 86-1(6-1), 86-2(11-6-1), 86-3(16-11-6-1), and 86-4(21- 
16-11-6-1) constitute the sequence of centered hexagonal packing. They are the most ideal 
cylindrical configurations in all the structures studied. In the double-chain-core sequence 
of wires 82-1(8-2), 82-2(13-8-2), 82-3(18-13-8-2), and 82-4(23-18-13-8-2), the central core 
consists a nearly parallel dimer-chain. The surrounding shells are made up of 8, 13, 18, and 
23 identical atomic strands respectively. Each strand helically wound with a linear pitch 
around the cylinder. The differences between the number of atomic strands in outer and 
inner shells are 5,5,5, and 6, respectively. The thickest wire in Fig.l is double-chain-core 
82-4 wire with D=2.787 nm. It is worthy noted that, in most growth sequence of zirconium 
nanowires, the number {n, n') of atomic stands that forming the inner and outer shells 
typically have even-odd coupling with n — n' = 5, with exception at 82-1 and 85-1 wires. 
These behaviors are similar to those obtained from previous experiment on suspended gold 
nanowire^. In the experiment, the number n, n' of atomic strands in inner and outer shells 
also exhibit even-odd coupling with n — n' = 1 , which can be related to the string tension. 
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Table II. Structures, diameters (D), average binding energies per atom (E), and average 
coordination numbers (CN) for the multi-shell zirconium nanowires with different growth 
pattern. The shell numbers n-nl-n2-n3-n4 for these wires is given in parenthesis. Bulk 
cohesive energy from present TB potential is Eb— -6.167 eV/atom. 
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The structures, diameters, coordination number, and binding energies of the zirconium 
nanowires in equihbrium structures are summarized in Table II. The average binding energies 
E and coordination numbers of zirconium nanowires generally increase as the diameter 
increase. The largest binding energy of the atoms in the thickest wires (S2-4) reach about 
95% of the bulk value (-6.167 eV/atom)i and that of the thinnest one (S3-1) is about 75%. 
If we plot the average binding energies E as function ofl/D, all the nanowires from different 
growth patterns fit well to the same linear dependence as: E = Eb — X/ D, where coefficient 
A=0.85 eV and the wire diameter D is in unit of nm. The linear dependence with 1/D can 
be understood by the surface effect0. The binding energies of all these multi-shell wires 
fitting to the same line indicate that they are almost equally stable. 
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FIG. 2. Angular correlation functions of zirconium nanowires. 

To illustrate the structural characters of zirconium nanowires, the angular correlation 
functions (ACFs) of these nanowires are calculated and presented in Fig.2. For most struc- 
tures, two major peaks in the ACFs located around 60° and 120° are found, which are 
related to the hexagonal packing. The broad distributions of ACFs in addition to the two 
major peaks demonstrate that noncrystalline multi-shell structures do not have any definite 
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bond angle. From the ACFs shown in Fig. 2, it is clearly seen that discrete feature of the 
ACF peaks gradually disappear as the wire diameter increases and the bulk character grad- 
ually come out. On the other hand, there are considerable differences between the angular 
correlation functions for the different growth sequences in Fig.2. For instance, the features 
related to deformed icosahedron can be observed in S5-x (x=l-4) sequence, while perfect 
icosahedral bond angles are 63.4°, 116.6°, and 18O°0. 
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Vibration frequncy (THz) 

FIG. 3. Vibrational densities of states of zirconium clusters. 

Based on the optimized structures of zirconium nanowires, we now their vibrational 
properties. The vibrational densities of states of zirconium nanowires in Fig. 3 exhibit re- 
markable dependence on the wire diameter and the corresponding growth classes. Except 
wires S3-1 and S4-1, the position and shape of the major vibrational peak located around 
2.1 — 2.7 THz do not sensitively change with the size and growth sequence of wires. This 
behavior is similar to that found for gold nanowires in our previous workEl. The first major 
peak around 2.1 — 2.7 THz can be related to the normal vibrational model of hep bulk 
zirconium^. In an earlier studyS, the frequency distribution of the bulk zirconium eval- 
uated at room temperature range between 1.2 — 6.1 THz. The first bulk frequency peak 
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were found at 2.1—2.3 THz, which is comparable to the major peak (2.1 — 2.7 THz) in 
zirconium nanowires. Different from the first major peak, although the positions of second 
peak (~ 5 THz) found for the nanowires do not substantially change, its shape and height 
sensitively vary with the size and growth sequence. Such differences might be understood by 
the different vibrational couplings between the atomic shells for the nanowires with various 
sizes and growth classes. In another words, the vibrational modes of zirconium nanowires 
sensitively depend on the specific growth patterns of the constituent atomic strands and 
diameters of nanowires. The vibrational bands of the thinner wires (S3-1, S4-1) are rather 
discrete and molecule-like. Some high frequency modes around l.bTHz are found in S2-x, 
S3-X and S4-x sequences. These high frequency vibration modes can be attributed to the 
contraction of interatomic distance in those wires, which will be further illustrated in the 
following discussion on zirconium clusters. 
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FIG. 4. Vibrational densities of states of zirconium clusters. 



In order to further analyze the relations between the vibrational spectra and atomic 
structures, we have calculated the vibrational densities of states for the zirconium clusters 
(Zry, Zri3, Zri5 , Zrig, Zr38, and Zr55) (Fig.4). These vibrational densities of states of clusters 
are comparable to those of the bulk and nanowires. Because of the structural similarity, the 
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shapes of vibrational densities of states for Zria and Zris are very closed to those of the 
pentagonal S5-1 and hexagonal S6-1 wires, respectively. The vibrational models of the 
clusters can be also understood from their atomic structures. The high frequencies are due 
to the large force constants, which are the consequence of the contraction of the interatomic 
distance in the clusters. For instance, the distances between surface atoms and central atom 
in an icosahedral Zri3 cluster (about 2.96—2.98 A) is smaller than these between surface 
atoms (3.10 — 3.14 A), which leads to the high frequency about 6.5 THz. Similar effect 
are observed in Zris, Zrig, Zrag, and Zrss. The appearance of some low- frequency modes for 
zirconium clusters may be attributed to the atoms weakly bonded to the cluster. 

From the above discussions, we can make the following conclusions. Multi-shell cylindri- 
cal structural growth sequences are found for the zirconium nanowires up to about 3 nm. 
The number of atomic strands in outer and inner shells exhibit even-odd and the differences 
between them is usually five for most growth sequence of zirconium nanowires. The angular 
correlation functions and vibrational properties are related to the wire diameter and growth 
pattern. The helical multi-shell metal nanowires represent a novel structural form of matter. 
They are similar but not identical to clusters, which is yet to be fully understood0ii. In 
principle, the mechanical, thermal, electronic, magnetic, and transport properties of such 
nanowires may sensitively depend on their atomic structures and helicity. The future stud- 
ies in these directions are under way. We expect future experiments on transition metal 
nanowires to validate our arguments. 

B. L. Wang and G. H. Wang would like to thank financial support from National Nature 
Science Foundation of China (No:29890210, 10023001). 
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